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' Abstract. In high energy nuclear collisions, QCD jets deposit a large fraction of their energy 

£SJ , into the produced matter. It has been proposed that as such matter behaves as a liquid with a 

very small viscosity, a fraction of this energy goes into a collective excitation called the "conical 
flow", similar e.g. to the sonic booms generated by the supersonic planes. In this work we 
study the effect of time-dependent speed of sound on the development of the conical wave. We 
show that the expansion of matter and the decrease of c s leads to an increase of observable 
manifestations of the conical flow. We also show that if the QCD phase transition is of the first 
order (and thus with vanishing speed of sound in the mixed phase) the wave must split into two, 
with opposite directions. We then argue that it is not the case experimentally, which supports 
the conclusion that the QCD phase transition is not of the first order. 
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Q_i! 1. Introduction 



A strong suppression of particle production at high transverse momenta has been observed at 
RHIC 1 , and it is commonly interpreted as due to parton energy loss (jet quenching). See early 
papers [2], more recent developments jSHUIS] and a short summary [H]. 
^ ■ An interesting question is where the deposited energy /momentum may eventually be 

observed, and in what form. Since the state produced at RHIC appears to be a liquid-like 
QGP [7j with a very short dissipative length (small viscosity) |5j, details such as the angular 
distribution of emitted gluons should be rapidly forgotten and energy /momentum be deposited 
locally. Further evolution would then be described by relativistic hydrodynamics. In 9 (in a 
brief form the idea was also mentioned in jlUj ) it was suggested that such local deposition of 
energy and momentum could start a "conical flow" of matter. The distortion of the shock waves 
due to the expansion has been studied in and an attempt to incorporate more realistic 
mediums has been performed in ,12 . Discussion of alternative explanations to the large pt 
correlations at RHIC can be found in ^3], [Tl], ^S] and [IB] . 

The solution of the linearized hydrodynamics equations (Hj provides a detailed picture of 
resulting flow of matter. The main direction of flow is at the so called Mach angle relative to 
the jet velocity 

cosOm = (1) 

Vjet 



where Vj e t is the jet velocity and c s is the (time averaged) speed of sound in the produced matter. 



For RHIC collisions the time-weighted c s till the freezeout time r is estimated to be 

cf HIC = 1/r f dtc s {t) « .33 (2) 

Jo 

Since for light quark and gluon jets Vj e t ~ c, the conical flow was calculated in |5] to be at about 
9 = 1.2 radian or at angles of about 70 degrees relative to the jet. Experimental observations 
made at RHIC, first by the STAR collaboration [Hj and then by PHENIX [IE], indicate a 
depletion of correlated particles in the direction of the quenched jet and a peak with an angular 
position and shape in agreement with hydro predictions. The issue whether this effect is real or 
follow from incorrect flow subtraction was discussed in detail by B.Cole ^B]) who demonstrated 
(slide 25) that the shape of the correlation function is quite independent of the direction of the 
jet with respect to the flow. Additionally one may observe it to be a real effect from the subset 
of the data, taken at the particular angle <j) = tt/4 at which the elliptic flow vanishes. These 
data show a clear minimum at Acf) = ir without subtraction. The maximum of the correlation 
function seems to be at an angle 110-120 degrees, an angular position in rough agreement with 
expected Mach cone position. The issue will be further clarified by the data on three-particle 
correlations reported at QM05: unfortunately at the moment their analysis is too uncertain to 
comment on them. 

In this paper we focus on phenomena which are induced by variable (time-dependent) speed 
of sound. This quantity, defined as 

c 2 = dp/de (3) 



via the thermodynamical variables pressure and energy density, is expected to change 
significantly during the process of heavy ion collisions. At early stages at RHIC the matter 
is believed to be in the form of quark-gluon plasma (QGP), and thus with Cq GP « 1/3. The 
next stage is the so called "mixed phase" in which the energy density is increasing much more 
rapidly than the pressure, so that c 2 decreases to the so called "softest point" and then 
rises again to c 2 RG = .2 in the hadronic "resonance gas". 

In this paper we will discuss basically two issues. The first is our finding that a gradual 
change of the sound speed affect the velocity amplitude of the wave. We apply rather well 
developed analytical methods based on adiabatic invariants (see textbooks e.g.{2Qj to infer how 
these changes 1 influence the final velocity of the sound wave, related the to final particles spectra 
in realistic conditions of heavy ion collisions. 

It is well known that for the description of the motion of any wave in a weakly inhomogeneous 
medium one can use geometric optics and eikonal equations to derive the phase of the wave. 
The amplitude magnitude is a more tricky question, which we would like to address with the 
help of adiabatic invariants 

I = <j> pdq (4) 

It is convenient to remind the reader what its conservation means for the basic example, 
the harmonic oscillator with a slowly variable frequency dlogio{t) / dlogt <C 1 and/or mass 
dlogM (t) I 'dlogt <C 1. The typical momentum and amplitude scale as 

p ~ yjM(t)E(t), q ~ ^E{t)/uj 2 (t)M(t) (5) 

Note also, that this analytic theory also helps to clarify a (somewhat methodical) question of monitoring 
accuracy of a numerical solution. In a stationary time-independent matter considered so far, this is done via the 
energy conservation, not available any more for any equations with time-dependent parameters such as c 2 . 



where E(t),tu(t), M(t) are slowly variable energy, frequency and mass. Their product is the 
adiabatic invariant I ~ E(t)/uj(t) which should remain a constant, independent on time. 

Hydro dynamical equations for sound waves, in momentum representation in coordinates, also 
form an oscillator, and therefore we will use the adiabatic invariant to study the changes on the 
velocity field due to the expansion and variable speed of sound. We will conclude that both of 
these phenomena enhance the sonic boom effect in QGP. 

The second issue is what will happen if the QCD phase transition is 1-st order, so that 
there is a truly mixed phase and the minimal value of c 2 is zero. As we will show below, 
this indeed would produce quite dramatic changes: all the waves, including conical ones, stop 
and then resume motion with splitting into two waves moving in opposite directions. Some 
phenomenological arguments at the end of the paper would lead us to believe that this in 
not what happens experimentally, in correlations with the high-pt jet. If this is confirmed by 
more accurate calculations and if the correlations observed at RHIC are due to hydro dynamical 
"conical flow" , it would mean that QCD phase transition can not be first order. 



2. Changing the conical wave amplitude by slow perturbations 

2.1. Sound in expanding matter 

Unfortunately, a simple substitution of a variable c s {t) into the equations of motion for 
perturbations of a static background is inconsistent. One should instead find a correct non 
static solution of the hydro dynamical equations and only then, using this solution as zeroth 
order, study first order perturbations such as sound propagation. The numerical solutions for 
hydro dynamical equations have been done by a number of authors: but in all of them the flow 
and matter properties depend on several variables and is too difficult to implement. 

Therefore, in order to study the effects of the variable speed of sound we have looked for 
the simplest example possible, in which there is a nontrivial time-dependent expansion but still 
no spatial coordinates are involved, keeping the problem homogeneous in space. The only way 
these goals can be achieved is by a Big-Bang gravitational process, in which the space is created 
dynamically by gravity. With such space available, the matter can cool and expand at all spatial 
points in the same way. For definiteness, consider a liquid in flat Robertson- Walker metric |21j : 

dr 2 = dt 2 - R(t) 2 \dr 2 + r 2 {d6 2 + sm 2 9#)] (6) 

where the parameter R(t) (the instantaneous Hubble radius of our "universe") is treated as 
external. This means that one should not consider Einstein equations of motion for it explicitly, 
just look at the effect of this gravity field on matter motion. 

The hydrodynamic equations (for simplicity, for a baryon free) fluid are the same as usual, 
but with the derivatives replaced by covariant derivatives (denoted by ;) 

T.»? = (7) 

The components of the stress tensor for ideal fluid are, as in the case of Minkowski space: 

= (e + p)u»u v - (8) 

This equation can be easily solved f° r a non-floating medium, = (1,0,0,0). It is 
straightforward to see that out of 4 equations the only nontrivial one is the longitudinal 
projection u^T. 1 ^ leading to the equation of entropy conservation 2 . Thus, 



dt 

li 

to the microscopic time scale. 



4 Ut)R 3 (t)) = =► s(t)R(tf = S (9) 



2 Note that this is also a variant of adiabatic invariant by itself, true if the matter expansion is slow in comparison 



As we already said, R(t) will be a conveniently chosen external parameter, that will drive the 
expansion dynamics. 

The next step is to derive the linearized equations for hydro dynamical waves in this 
background metric 8Tf^ = that we can express as, 

R 



d t ^(-g)5T 00 j + d t (y (-0)*Z*J + 3-c^T 00 v /(- 5 ) = (10) 



Where we have use c s = dp/de in order to relate 5T lJ to ST , and \f\—g) = Br. We define 



e = R 4 5T 00 (11) 
(? = R 5 5T 0i 

and perform the standard change of variables 

V = J^dt (12) 

such that the metric takes a conformally flat form. After some straight forward algebra one gets 
the equations of motion for small perturbations: 

d v e + diG l + {3c 2 s - l)^e = (13) 
R 

drjC? + €%€ = 

where here R' = dR/drj. 

In the derivation of (jl.3[) we have not assumed any particular expansion law or equation of 
state, just thermodynamic properties. Given the time dependence of the Hubble radius s(t) and 
provided the equation of state for the fluid, we can obtain all the thermodynamic properties, in 
particular c s . 

In order to get a better feeling about the different contributions in (|13|) let us write the second 
order equation for the quantity e 

d 2 v e - c 2 s Ve + ed v ((3c 2 s - 1)^) + (3c 2 - l)^d v e = 0. (14) 

The resulting wave equation has two new terms, the first changes the speed of wave propagation 
and the second is a rescaling factor. Note that the corrections to the wave equations vanish for 
c 2 = 1/3 (the QGP value) and the rescaling factor actually produces an amplitude growth, if 
the speed of sound reduces below the ideal value c 2 < 1/3. 

Let us Fourier transform the space components and introduce the following quantities: 

Tr = ((3^-i)-n) ( 15 ) 



M V R 



LO' 



d v (3c; 



4 



(16) 



then, we can write (|14j) as the equation for a harmonic oscillator with (time dependent) mass 
M and frequency lo 



d r] (Md rj e)+Moj 2 e = (17) 
As is well known, the Hamiltonian for such system can be written in the oscillator form 

with p = Mde / dn and q = e 

One can do a similar analysis for the momentum G. As in in order to have propagating 
modes (sound waves), the space Fourier transform of G should be in the form: 

G(t,k) =kG L = ikG (19) 

With this requirement we can also write a harmonic oscillator kind of equation for G that looks 
slightly different from the previous. After Fourier transform of 1)13(1 the second order equation 
for G is 

d v (Md n G) + Mlu 2 G = (20) 

with M = M/c 2 and Q = kc s . The Hamiltonian H as well as the generalized coordinates q, p 
are defined in an analogous way to the previous case. 

2.2. Adiabatic invariants 

Adiabatic invariants allow to predict properties of open mechanical systems subject to small 
(adiabatic changes) in the external parameters. We will now use those to learn about the effect 
of the expansion in the final solutions as well as to check our numeric procedures. After the 
analogy of the equations for the (space) Fourier modes with those of the harmonic oscillator 
with time dependent mass and frequency is established, it is rather straightforward. 

The adiabatic condition for any time dependent parameter A of the equation of motion (M 
or u) in our case) is 

where T osc is the period of oscillation of the solution of the equation of motion of the system 
considered as if the parameters where constant. In our case T osc = 2ir/u 

As stated before, if the parameters change in an adiabatic fashion, the adiabatic invariant 
/ remains approximately constant through out the evolution of the system. For the harmonic 
oscillator the following identity holds 

/ = ^ (22) 

u> 

where H is the Hamiltonian and the average is taken though one period of oscillation of q (taken 
as if both M and u> where constant and equal to the value at some time). Thus, the amplitude 
Ak for an oscillation e(k, t) is 



(23) 



In the same way, one can perform the same analysis for the time dependent harmonic oscillator 
for G and obtain for the amplitude of the oscillations 



Note, however, that these quantities are related to the components of the perturbed stress 
energy tensor though time dependent quantities. In particular we would like to look at the 
behavior of the perturbed velocity field v l = Ru % (where the factor R recovers the correct 
dimensions), related to 5T° l = wu 1 , with w the enthalpy w = e + p = Ts. Thus, we find 

(~t% (~ti 

' (25) 



R A w STR 

We now can relate TR to M remembering that for a baryon free fluid c 2 s = s/TdT/ds, and by 
use of equation © 

V 9 R' ( M' R'\ , . 



(27) 



T 6 R \M R 

Thus, we can fixed the normalization for M such that 

1 



MR 

Combining this with the previous expression, 

,■ G'M 

v = 

S 

And thus the amplitude of the velocity (v\ for a fixed Fourier mode ) behaves as 



(28) 



VM02I (29) 



Therefore, for a medium with c 2 s < 1/3, where M decreases with increasing R ()15|) . we conclude 
that the amplitude of the perturbations of the velocity field decreases with the expansion. The 
dependence on the speed of sound is two fold, through Q and M. From (|15[) we see that for a 
fixed expansion rate and initial value for the temperature, M decreases faster for smaller c s . As 
also u> decreases in an obvious fashion with c s , the amplitude of the velocity perturbations also 
decreases with decreasing speed of sound. 

However, as we will show when we discuss the spectrum, the relevant quantity for the final 
production of particles is not the velocity but the ratio v l /T = Rv l /M. Thus, the ratio of the 
amplitude of the wave (with mode k) to the temperature behaves as 



as in the previous case, for (? s < 1/3, it increases with the expansion. The dependence on 
the speed of sound is more complicated due to the ratio u>/M, where both quantities decrease 
with c s . So in order to determine the dependence on c s we need extra information on the time 
dependence of both quantities. 



10 

10* 



10 

10 ' 
10' 
10' 
10"" 



i 0.16 ; 
< 0.14 ; 

o.i2 ; 
o.i ; 

o.oa ; 

o.o6 ; 

0.04 ; 

0.02 : 



80 100 120 140 160 180 



Figure 1. Check of our numerical procedure for the simplest possible case (see text) (a) 
time dependent mass M for the k=2 mode (b) Relative change of the adiabatic invariant 
AI/I = (I(r]) — I(0))/I(0) in 20 periods of oscillation. In both cases 7 = 0.1, 0.5, 1 for 
solid, dashed and dotted lines respectively. 



2.3. Monitoring the numerical solution 

In the next section we will study a case where the adiabatic approximation is not valid. Then, 
we will need to solve numerically the system of equations ((T3|) . We use the point splitting 
technique, separating (|13j) in a wave equation and a rescaling part. For the wave equation we 
use the MacCormack's technique |22j . and for the rescaling factor a first order Runge-Kutta 
update. 

The adiabatic invariant provides us with a way to check our numerical procedure. In the 
static medium case, one can use the energy conservation to monitor the solution. Now that 
check is not available and, thus, we need to look for some other (approximate) constant of 
motion, that is the adiabatic invariant. 

We will illustrate our check in the simplest possible case (after c 2 s = 1/3), that of constant c s 
and constant frequency of oscillation mode by setting 

R = e ari (31) 
with arbitrary a . With this prescription, only M depends on 77 in the following way 

M = d 3 ^- 1 ) (32) 
The adiabatic condition for this system reads 

M' u 3c? - 1 

— « — -> 7 = 2vr s , a « 1 33 
M 2tt c s k y ' 

In Fig^we plot the relative change of the adiabatic invariant (AI/I = (I(rf) — 1(0))/ 1(0)) for 
a particular Fourier mode (k=2) for different values for the coefficient 7. In practice, we fixed 
the value of 3c^ = 0.5 and changed a. In Fig^a we show, the change of M for different choice 
of parameters. In Fig E]b we show the change of I. Even though the mass changes by orders 
of magnitude, our numerically calculated adiabatic invariant remains constant on the level of 
0.01 %. Even for the values of 7 = 1 the total change of the adiabatic invariant in the time 
considered (20 T osc ) is very small. 



2-4- The wave amplitude and the observed spectrum 

The spectrum of particles produced by hydrodynamical velocity and pressure fields is determined 
at the freeze out. The standard way to calculate it is the Cooper- Fry prescription, and we discuss 



how it works in the case we discussed, with cosmological-style expansion. Although it is not 
identical to experimental conditions, with flow fields depending on spatial coordinates, one can 
still get some lesson from this example. 

The simplest way to obtain a particle spectrum from the hydro calculation is to consider 
that the system behaves hydrodynamically till some freeze out happens and after that particles 
(created through a thermal distribution) free stream. If the case of flat space, when the spectrum 
is calculated though the Cooper-Fry prescription at some time, the free streaming stage does 
not change the invariant spectrum. 

However this is not the case in our curved expanding space because particle momenta 
continue to change in the expanding Robertson Walker space (as an illustration, the cosmic 
microwave background has a Boltzmann distribution, but the temperature of the distribution 
is time dependent). Thus, if we would keep the radius of our "universe" to expand for ever, 
the spectrum of particles we calculate after freeze out will continuously change, dropping the 
measured temperature. 

In order to obtain an answer that resembles more the experimental conditions one may 
consider the "universe" which stops expanding after the freezeout. Thus, the spectrum remains 
constant when the particles fly out of the collision region. As in [S] we will consider fixed proper 
time freezeout, for which 



,„ = I R 3 f d 3 x 1, (34) 

Where Rf is the value of the universe size at freeze out (and constant afterwords as an 
assumption). This factor has to be included to get the correct dimensions. The exponent of the 
thermal distribution also includes such factor. For a small perturbation, it will be 

pPup = E - R^RfU 1 (35) 

From here we observe that for a fixed momentum (Rtp 1 ) the relevant quantity that sets the 
correlations in the spectrum is v/T. We have shown already that this combination of fields 
increases with time. Thus, the effect of a small disturbance at early times, when the medium is 
more dense, increase as the system expands and becomes more dilute, as the amplitude of the 
relevant hydrodynamic field grows. 



3. Vanishing speed of sound and reflected waves 

As stated in the introduction, the speed of sound reaches a minimum in the mixed phase. If 
the system experiences the first order phase transition, the speed of sound in the mixed phase 
should vanish, because the pressure does not change during this stage while the energy density 
does. In this section we study the consequences of a vanishing speed of sound for the sound 
propagation. 

Let us start by noting that when the speed of sound vanishes, the adiabatic approximation 
explored in the previous section can no longer be applied. In fact, in such a case, the logarithmic 
derivatives of c s diverges and thus, (|21|) does not hold. Thus, we will need to solve numerically 
the equations ()13|) with the procedure checked in the previous section in order to investigate 
this effect. 

In order to get somewhat realistic estimates of the effects, we will model both the expansion 
rate and the speed of sound according to the physical situation at RHIC. Inspired in the boost 
invariant scenario, we will chose R(t) such that the entropy density behaves as: 

s(t) = s 0J ^- t R(t) = R (^) V3 (36) 
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Figure 2. (a) c s as a function of t. The two curves are different interpolations with different 
minimal at the mixed phase (0 (solid) and l/2\/2 (dashed)) (b) Ratio of amplitudes of the 
reflected and transmitted waves as a function of the minimal c s in units of c s of the hadron gas 

Where to = 0.5/m and Rq = 6fm. The value of sq will be specified later. 

We now discuss the function c s (t) . Numerical solutions of hydrodynamics as in |23| use values 
of the speed of sound that are discontinuous at the beginning of the three phases. Here, however, 
we study parametrization of the speed of sound that interpolate smoothly between the three 
stages of matter at RHIC. The times for the interpolation are also inspired in hydrodynamic 
solutions for RHIC j2B] • Following these calculations we will assume that each of the three stages 
have the same time duration of 4 fm. Finally, in order to investigate the effect of the first order 
phase transition on the sound waves, we will leave the minimum value of the speed of sound in 
the mixed phase as a free parameter and we study the propagation of the waves for different 
values of this minimum, (see Fig. Ela) )• Once the function c s (t) is determined, we can solve 
(|2fij) to obtain the temperature, imposing that at the time where c s (t) has a minimum (t m i n = 6 
fm in our case) T(t m ^ n ) — T c . This determines uniquely the temperature, and in particular the 
value T(0). We use this value of the temperature and the equation of state for an ideal QGP in 
order to determine so- 

Once all the parameters are specified, we proceed to numerically solve (j!3j) . We will study 
the propagation of a spherical wave started from a perturbation at time t=0. Let us start our 
discussion with the case of vanishing minimum c s (first order phase transition). In this case, the 
spherical wave will split into two at the mixed phase, one propagating outwards (transmitted 
wave) and the second inwards (reflected) towards the origin. A simple physical explanation 
of this effect is as follows: when the speed of sound vanishes, all the information about the 
directionality of the wave disappears and, thus, when the speed of sound changes again to a 
finite value the disturbance splits into two, as it can propagate both inwards and outwards. 
Finally, the spherical wave propagating inwards eventually bounces back when it arrives at 
the origin, leading to a second outward moving spherical wave. This is observed in numerical 
solution as a double peaked structure of the spherical wave profile at sufficient late times. 

The appearance of the reflected wave happens not only when the speed of sound vanishes. 
We studied the spherical wave profile for different values of the minimum. This is best done 
at sufficiently late time, when the two waves (transmitted and reflected) propagate in the same 
direction. To illustrate the importance of the reflected wave, we study the amplitudes of the 
maximums in the wave profile for both waves at a late time. (In order to get rid of a trivial \/r 
decay of the wave, we multiply them by the position of the maximum r.) Thus we define the 
relative amplitude: 



T ^reflected 

f '^transmitted 



(37) 



where the values are to be taken at the maxima of the disturbances. We plot this quantity in 
Fig El b). We observe that, as claimed, the reflected wave appears for non vanishing values of 
the minimal speed of sound. However the relative importance of the reflected wave decreases as 
the minimal value increases (as expected) 

Even though the particular dependence of the relative amplitude with the minimal speed of 
sound, depends on our assumed parametrization c s (t), the general phenomenon is independent 
of the parametrization. As long as the speed of sounds vanishes at the mixed phase, we will 
obtain a reflected wave for the spherical disturbance. 

Finally, for the cases where the minimal speed of sound does not vanish (as in the dashed 
curved in FigEJa ), we can use the adiabatic invariant to learn about the amplitude dependence 
of the disturbance. As we have specified a particular time dependence of R(t) and c s (t), we can 
complete the analysis of the previous section for this particular case. By solving (|15|) we find that 
for our parametrization of the speed of sound with minimum ranging from 0.5\/0.2 to 0.9\/0.2, 
the effect of the overall dropping speed of sound increases the ratio v/T by 10%. As before, the 
expansion of the medium also induces an amplitude increase. For our particular expansion rate, 
setting the freeze out time to tt ~ 14/m, we obtain for the same quantity a factor of 3 increase 
(as the system becomes three times larger). These changes are very important experimentally, 
as the ratio v/T enters in the exponent of the Cooper-Fry formula for the observed spectra. 

3. 1 . Comparison to experimental correlation functions 

Now we discuss where such reflected waves, if produced, will be located in the conical flow 
induced by high energy jets. The cone arises from the superposition of spherical waves produced 
along the path of the jet in the medium. All those disturbances produced prior to the mixed 
phase will encounter a time of minimum speed of sound. Following our previous discussion, if 
this minimum is small enough, all those waves will split into transmitted and reflected waves. 
After the superposition, given sufficient time, the reflected waves will generate a second cone. 
So, provided the observation time is late enough as for all the reflected wave to bounce back, 
one would find two cones with different opening angle propagating in the same direction. 

However if the observation time is too close to emission time, one would find the the second 
cone moving inwards (or opposite direction to the jet). At RHIC the timing of the expansion 
is such that the three main stages - the QGP, the mixed phase and the hadronic phase - take 
roughly comparable period of proper time, about 4-5 fm/c each. The speed of sound at the 
QGP phase c® GP = l/y/3 is larger than that in the final hadronic phase = \J~H. Therefore, 
the distance the conical wave would propagate in QGP is larger than that in hadronic phase, 
and thus the reflected wave cannot possibly reach the origin (for spherical wave, the original 
axes for conical one). 

What this conclusion implies experimentally, is that if the QCD phase transition is first 
order and the original wave gets split into direct and reflected waves, the freezeout will find the 
reflected wave on its way to the origin, the opposite to the Mach direction of the direct wave. 
In terms of two particle azimuthal distribution it means the appearance of a second peak for 
A(p < it/2. Even though detailed numerical simulations are needed to give the correct position 
of the peak, we give here a simple estimate for its location, based on geometrical optics. 

The disturbance generated at the origin travels a distance AB 

ftM 

AB = / dtc s (t) (38) 
Jo 

before the minimum of c s is reached and the wave splits into two. During the mixed phase 
the disturbance does not propagate, but after the hadronic stage starts, the waves propagates 




Figure 3. Azimuthal dihadron distributions normalized per trigger particle measured by 
PHENIX [IH] for two different centralities (2.5 < p ^ er < 4GeV,l < p f sociated < 2.5GeV). 
The filled arrow indicates the position of the Mach cone. The empty arrow our estimate for the 
position where the cone from the reflected wave should appear. 



backwards a distance 

BC = f 1 dtc s {t) (39) 

Where tf is the freeze out time and we have assumed that the mixed QGP phase and the 
mixed phase have the same time duration. Thus, the relevant distance traveled by the wave is 
AC=AB-BC. The relevant longitudinal distance is not the distance traveled by the jet, but the 
distance traveled till the minimum of c s , ijv/ (as after that time there is not more "emission" of 
reflected waves). Thus, the effective angle of emission is 

cos(9 e ) = AC jt m => e ~ 1-4 rad (40) 

where we used for AB and BC fixed c s values of 1/V3 and \/0.2 respectively. 

In Fig. EH we show samples of experimental correlation function as obtained by PHENIX |18j . 
The peak around 1.9 rad (indicated by the filled arrow) corresponds to the Mach direction; this 
is the main effect attributed to a conical flow. The reflected wave should appear in the region 
A(p < 7r/2, and the empty arrow shows the estimated place where the corresponding reflected 
peak should appear. Fig|3]a) shows the most central collisions, does not show any nonzero signal 
at that angle. Fig|3]b), at higher centrality, has a nonzero correlation function there, but it 
seems likely to be just a slope of a much broader peak. 

We thus argue that there are no indications for enhanced correlations at the expected angle 
of the reflected wave, and thus the deconfining phase transition cannot be of the first order. 

4. Summary and discussion 

In this paper we have studied the effect of a variable speed of sound in the propagation of sound 
waves. We have considered for simplicity an expanding Universe with a flat Robertson- Walker 
metric, that allows us to introduce dynamical expansion depending on proper time only. 

We have use the adiabatic approximation to investigate the effect of the expansion on the 
waves, concluding that even though the velocity field decreases with the expansion, the ratio 
of the velocity to the temperature increases. This ratio enters the exponent of the Cooper-Fry 
formula, and thus the observed effects in particle spectra may be significantly enhanced, the 
more so the larger is the particle momenta. We have shown that within the model studied based 
on the Robertson- Walker metric we obtain a 10 % increase of the ratio v/T due to the change 



of c s and a factor 3 increase of this ratio due to the expansion. Even though this values depend 
on the particular model studied we expect a similar effect to take place in fireball expansion. 

We have also studied the effect of a vanishing speed of sound, where the adiabatic 
approximation cannot be used. This is supposed to happen in the mixed phase if the phase 
transition in QCD is first order. We have shown that this must lead to the appearance of 
reflected waves, the superposition of which leads to a second cone in the flow field. Given the 
timing of the different stages at RHIC the second cone should be moving inwards at the time of 
freezeout. We concluded that the first order phase transition should lead to a peak at a relative 
angle Acfi rs 1.4 rad. The experimental correlation functions do not show any enhancement at 
this direction, which we think indicates that the phase transition in QCD cannot be of the first 
order. Needless to say, more theoretical and experimental work is needed before finalizing this 
important conclusion. 
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